I-V curves of Fe/MgO (001) single- and double-barrier tunnel junctions 
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In this work, we calculate with ab initio methods the current-voltage characteristics for ideal 
single- and double-barrier Fe/MgO (001) magnetic tunnel junctions. The current is calculated in the 
phase-coherent limit by using the recently developed Smeagol code combining the non-equilibrium 
Green's functions formalism with density functional theory. In general we find that double-barrier 
junctions display a larger magnetoresistance, which decays with bias at a slower pace than their 
single-barrier counterparts. This is explained in terms of enhanced spin-filtering from the middle Fe 
layer sandwiched in between the two MgO barriers. In addition, for double-barrier tunnel junctions 
we find a well defined peak in the magnetoresistance at a voltage of V = 0.1 Volt. This is the 
signature of resonant tunneling across a majority quantum well state. Our findings are discussed in 
relation to recent experiments. 

PACS numbers: 85.75.-d, 72.25.Mk, 73.40.Rw, 73.23.Ad 



I. INTRODUCTION 



Tunneling magnetoresistance (TMR) arises in mag- 
netic tunnel junctions consisting of an insulating barrier 
sandwiched by two ferromagnetic electrodes. It is ob- 
served that the resistance strongly depends on the rela- 
tive orientation of the magnetization vectors of the elec- 
trodes. It is usually small when the magnetizations are 
parallel to each other and increases when they are an- 
tiparallel, and the change in resistance can be large. The 
TMR coefficient is defined as TMR = [(I P -I AP )/I AP } x 
100, where I P and Jap are the currents in the parallel 
(P) and in the antiparallcl (AP) magnetic configuration 
related to the same applied voltage, respectively. Since 
the pioneering work of Julliere [J HI, the experimentally 
attainable TMR values have been steadily increasing, 
mainly due to the tremendous advances in growth tech- 
niques. In particular a steep increase in TMR magnitude 
followed from the growth of highly crystalline Fe/MgO 
(001) tunnel junctions with atomically flat surfaces Q. 
Currently the record value is 500 % at room tempera- 
ture in fully epitaxial FeCoB/MgO/FeCoB junctions [4(. 
This value is much larger than the one it can be obtained 
in standard metallic giant magnetoresistive devices, and 
this explains why magnetic tunnel junctions will soon re- 
place metallic spin- valves in the read heads of near-future 
hard disk drives ||. 

In epitaxial single-barrier magnetic tunnel junctions 
(SBMTJs) the large TMR is a consequence of the sym- 
metry matching between Bloch states in the electrodes 
and the complex bands within the insulator [6]. These 
latter are evanescent states lying in the band gap and 
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decaying inside the spacer. For thin spacers, another 
transport mechanism may arise, when surface states at 
each side of the barrier come into resonance. These form 
the so-called hot spots in the transmission coefficient as 
a function of the transverse wave-vector fen. This mech- 
anism was clearly demonstrated in the ab initio calcu- 
lations of MacLaren et al. 7] for Fe/MgO SBMTJs, in 
which a Fe minority surface state was found resonant 
through the MgO spacer. It was then shown experimen- 
tally by Tiusan et al. [8| and theoretically by Rungger 
et al. 9j that, in Fe/MgO SBMTJs, these surface states 
are very sensitive to the applied bias voltage V, and as a 
consequence the TMR is maximal at zero bias and signif- 
icantly decreases with increasing bias. The bias voltage 
at which the TMR drops to half of its zero bias value 
is denoted as Vi/% and it is taken as one of the quality 
factors for the applicability of SBMTJs in real devices. 
Typically Vi/ 2 is less than 0.7 Volt. 

Nowadays, one of the main challenges in the field of 
spin-electronics is to reach simultaneously large TMR 
and Vi/2 values. One possible route to accomplish this 
goal was proposed by Zhang [llj a decade ago and con- 
sists in the insertion of a metallic slab in between the 
insulating spacer, thus to form a double-barrier magnetic 
tunnel junction (DBMTJ). In this architecture, quantum 
well states (QWSs) formed by confinement are found in 
the in-between metallic slab. These can come into reso- 
nance with the evanescent states inside the barriers, en- 
hancing the TMR in a significant way and allowing the 
tuning of the TMR by controlling the thickness of the in- 
between slab. This tunability of the TMR arises because 
changing the thickness of the in-between metallic slab 
produces a shift in the energy of the QWS, and creates 
new ones. In addition, since the potential drop is shared 
across the two barriers and an additional spin-dependent 
scattering potential is introduced by the in-between mag- 
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netic slab, the V\ii for DBMTJs is expected to be larger 
than that of a SBMTJ of same combined thickness. 

To date there have been several theoretical studies on 
the spin-dependent properties of DBMTJs p], [H [H, 
Hij . with either magnetic or non- magnetic in-between 
slabs. Usually these are conducted at zero bias in the 
linear response limit or in a non-self-consistent fashion. 
The key features emerging from these studies are, first 
that the TMR can reach extremely large values under 
resonant conditions, and secondly that the TMR can be 
enhanced not only by resonant tunneling through QWS 
but also by the spin-filter effect (SFE) Q. This effect 
is a consequence of the insertion of a magnetic slab in 
between the barriers. In fact in these DBMTJs the P 
configuration corresponds to all the magnetizations (the 
two electrodes and the in-between layer) being aligned 
parallel to each other, while in the AP configuration the 
electrode magnetizations remain parallel to each other 
but the in-between magnetization is antiparallel to them. 
Thus in the AP configuration of DBMTJs there are ef- 
fectively two interfaces where the magnetization vector 
change sign and this greatly enhances backscattering. 
Since the electrodes and the in-between Fe layers have 
different coercive fields, due to their different thicknesses, 
these magnetic configurations P and AP are experimen- 
tally attainable [f|. 

The SFE has been clearly shown in our previous cal- 
culations for ideal Fe/ZnSe_{001) SBMTJs and DBMTJs 
with Fe in-between layers [14j- The main difference be- 
tween resonant tunneling through QWS and the SFE is 
that the former strongly depends on the thickness of the 
in-between slab, while the latter is almost independent. 
Thus, DBMTJs show a variety of spin-dependent trans- 
port phenomena richer than that of their single-barrier 
counterparts, and in principle they promise a better con- 
trol of the TMR. Some of these expectations are now 
confirmed experimentally. For example, Nozaki et al. 
[lH have recently measured the I-V curves of fully epi- 
taxial Fe/MgO SBMTJs and DBMTJs (with Fe as the 
in-between metallic slab), and found that (i) the TMR 
of DBMTJs is larger than that of SBMTJs, and (ii) 
the TMR decrease with bias is significantly slower in 
DBMTJs than in SBMTJs. Similarly Zeng et al. [HI 
have shown that CoFeB/Al-0 DBMTJs also have TMR 
and V1/2 values larger than those of the corresponding 
SBMTJs. Other works along the same lines which con- 
firm these results include that of reference [ItJ ■ 

In spite of the amount of research carried out on the 
spin-dependent transport properties of DBMTJs, self- 
consistent calculations of I-V curves for a realistic junc- 
tion from first principles are still lacking. The purpose 
of this paper is to fill this gap. We present ab initio cal- 
culated I-V characteristics for Fe/MgO (001) DBMTJs, 
and compare them to those of SBMTJs, relating the 
transport properties of the devices to their electronic 
structure. To this goal, we use the recently developed ab 
initio code Smeagol 18], that combines the pseudopoten- 
tial density functional code Siesta [19( for the electronic 



structure calculations with the non-equilibrium Green's 
function formalism (NEGF) for phase-coherent transport 

sa- 
il. CALCULATION DETAILS 

Our SBTJs consist of n monolayers (MLs) of MgO 
(001) sandwiched by two semi-infinite bec Fe (001) elec- 
trodes, while our DBMTJs are multilayers of the type 
(MgO)„/Fe m /(MgO)„ (001) sandwiched by the same 
electrodes. In both cases, the junctions are assumed to 
be periodic in the x-y plane, being z the transport di- 
rection. In order to account for the charge transfer and 
to correctly reproduce the band offset between Fe and 
MgO, we include in the cell for self-consistent calcula- 
tions four Fe MLs belonging to the electrodes at both 
sides of the junction. This is enough to correctly account 
for charge screening inside the ferromagnet. Similar to 
previous calculations @, H, HH the lattice constant of the 
electrodes is fixed to 2.87 A and that of MgO is taken to 
be v2 larger. This, together with a 45° rotation of the Fe 
unit cell, allows epitaxial matching between Fe and MgO. 
Figure [T] shows the structure of a single-barrier junction 
with n=2 MLs (4 A), together with the schematic struc- 
ture of single- and double-barrier junctions. In this work, 
the possible appearance of FeO interfacial layers, as well 
as atomic relaxation and disorder, are not considered, so 
that our calculations are valid only in the ballistic limit 
for atomically ordered structures. 

For the electronic structure of the junctions, we use 
norm-conserving pseudopotentials, double-zeta basis set 
for all the angular momenta and the generalized gra- 
dient approximation (GGA) [2l[ to the exchange and 
correlation potential. We have thoroughly checked that 
the band structure and the density of states of bulk Fe, 
bulk MgO and Fe/MgO multilayers, as well as the charge 
transfer and magnetic moments in the last case, are very 
well reproduced as compared to FP-LAPW results ob- 
tained using the highly accurate WIEN2k code [22[ ■ We 
obtain a band offset (the difference between the Fermi en- 
ergy E-p of Fe and the valence band of MgO) of 3.51 eV, 
in very good agreement with previous theoretical [23| and 
experimental reports [24| . As well-known, density func- 
tional calculations using semi-local exchange and corre- 
lation functionals underestimate the band gap and ours 
are not an exception. We obtain a band gap of 5.4 eV (as 
compared to the experimental value of 7.8 eV [24 J . 12511) . 
which agrees well with what expected from GGA p], [23| . 

The ballistic current density at each bias voltage V is 
calculated as 

I a (V) = I J dE T°{E, V)[f L - /r] (1) 

where / L = f(E - fj, h ) (f R = f(E - /x R )) is the Fermi- 
Dirac function evaluated at E — /Zl (E — /x R ) with fij, = 
E F + eV/2 (/i R — E F — eV/2) the chemical potential 
of the left (right) electrode. Finally a is the spin index 
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FIG. 1: (a) Structure of single- and double-barrier junctions. 
The arrows indicate the two magnetic configurations consid- 
ered (parallel/antiparallel). (b) Structure of a SBMTJ with 
n=2 MLs, showing the active region and part of the semi- 
infinite Fe electrodes. The junction is periodic in the x-y 
plane while z is the transport direction. Distances are in A. 



standing for majority (a =|) and minority (a =J,) spins. 
The transmission coefficient T cr (E,V) is calculated for 
each bias and it is given by 



T°{E,V) = [ dk x dk y T°(E,V,k x ,h y ) 
vbz J 



(2) 



where Vbz is the area of the 2D Brillouin zone orthog- 
onal to the transport direction. Here we assume that 
both spin and transverse momentum are conserved, an 
approximation that is valid for relatively thin epitaxial 
junctions. The k^-resolved transmission coefficient ap- 
pearing in Eq. (fc|| = k x x + k y y, see Fig. [1]) is 
calculated from the non-equilibrium Green's functions 
formalism in the usual way [l8|, |2(|. It is given by 
T = Tr[Ti,G T T^G a }, where for simplicity we omit the 
spin label a. Here, Ti^r are the broadening matrices de- 
scribing the interaction (thus the finite lifetime) of the 
scattering region energy levels with the left- and right- 
hand side electrodes, and G r (G a ) is the associated re- 
tarded (advanced) Green's function describing the one- 
electron electronic structure of the scattering region. The 
broadening matrices are calculated from the self-energies 
^l,r a s Tl,r = i{^L,R — R ). These in turn are ob- 



tained with the semi-analytic method described in refer- 
ence dHH]. 

In our calculations, we use a 8 x 8 x 8 k-point mesh 
in reciprocal space to calculate the density matrix of the 
scattering region and a 150 x 150 x 1 mesh to evaluate the 
current at each bias voltage. We have carefully verified 
that these meshes are sufficient for converging the density 
matrix and the current. 



III. RESULTS: SINGLE-BARRIER TUNNEL 
JUNCTIONS 

In order to benchmark our calculations for DBMTJs, 
we have calculated first the conductance and TMR at 
zero and finite bias of several single-barrier Fe/MgO (001) 
junctions as a function of the barrier thickness n. In 
general we find a rapid decrease of the TMR as a function 
of bias in agreement with previous theoretical results [j| . 
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FIG. 2: Current density and TMR as a function of bias volt- 
age of the single-barrier junction with n—2 ML. 

As an example Figure [5] shows the current-voltage 
curve for n=2 ML (4 A), together with the correspond- 
ing TMR. At zero bias, the TMR is calculated from the 
transmission probabilities, Eq. (J2), evaluated at E-p. It 
is seen that the currents for each spin channel are al- 
most linear functions of V over the bias range inves- 
tigated. However, the P minority current saturates at 
V = 0.4 Volt while the P majority and the AP currents 
keep increasing beyond that value. The competition be- 
tween the rates of increase as a function of bias of the 
P majority and AP currents produces a rapid decay of 
the TMR and we calculate V1/2 being at around 0.7 Volt. 
The TMR peaks at zero bias where it reaches up to 290 
% but, because of its almost exponential decay, it ap- 
proaches the mean value over the bias range investigated 
of 160 % within only 0.2 Volt. All of these features are 
in good agreement with experimental data [l!, [H, [29[ 
and early theoretical results [3(| HH . It is important to 
note also that, at very low bias voltages, the P minority 
current is slightly larger than the AP current. As it was 
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shown by Rungger et al. 0] , this large P minority current 
near to zero-bias is due to the resonance of a Fe minority 
surface state through the MgO barrier and it is washed 
out as soon as a voltage is applied across the device. 



IV. RESULTS: DOUBLE-BARRIER TUNNEL 
JUNCTIONS 

Figures [3] and 3] show the current density and the TMR 
for junctions with n=2 MLs and m=2 MLs (2.87 A) and 
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FIG. 3: Current density and TMR as a function of bias volt- 
age for a double-barrier magnetic tunnel junction with n—2 
MLs and m=2 MLs. 



n = 2, m = 4 
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FIG. 4: Current density and TMR as a function of bias volt- 
age for a double-barrier magnetic tunnel junction with n—2 
MLs and m=4 MLs. 

tice is the appearance of oscillations in the TMR at low 
bias. These are not present in the SBMTJ, and stem from 
peaks in the low bias I-V curve for all the spin channels. 
It is seen that the P majority current has a large peak 
at V — 0.1 Volt, present for both values of m. Although 
the P minority and the AP currents also present peaks 
as a function of bias, these are considerably smaller than 
those in the P majority channel. A difference between the 



two junctions investigated is that the bias at which the 
peaks appear depends on the number of monolayers m of 
the in-between Fe layer. Thus, the TMR oscillations as a 
function of bias are different in the two cases. In particu- 
lar, for m=2 there are two TMR peaks at V = 0.05 Volt 
and at V = 0.15 Volt, at which the TMR reaches values 
close to 1250 %. At V = 0.1 Volt, the P majority current 
peak competes with that of the AP configuration, result- 
ing in the actual decrease of the TMR. In contrast, for 
m=4, there is a very large TMR peak at V = 0.1 Volt, 
that clearly originates from the corresponding one in the 
P majority current. 

The peak in the P majority current of DBMTJs is 
suggestive of resonant tunneling through a QWS posi- 
tioned in the in-between Fe slab. This is consistent with 
the experimental reports of Nozaki and coworkers [32| 
for Fe/MgO (001) DBMTJs with Fe islands as the in- 
between slab. They observed small conductance oscilla- 
tions with bias only in the P configuration. In particular, 
for the thinnest in-between Fe slab investigated (1 nm), 
they observed two peaks in the P differential conductance 
at V ~0.15 Volt and V -0.4 Volt. Very recently, Wang et 
al [HI have theoretically shown that these conductance 
oscillations could be fit to majority QWSs of Ai sym- 
metry at the r point, once the level shift produced by 
charging effects is taken into account. As it can be seen 
in Fig. 2(a) of reference for ultrathin in-between Fe 
slabs there is a majority QWS at the T point with an en- 
ergy slightly less than Ep+0.1 eV. It is then tempting to 
ascribe our P majority peak at V =0.1 Volt to resonant 
tunneling through the above mentioned QWS. 

Although our intention in this work is not to push for 
a quantitative agreement with the results of Nozaki et al 
(321 ] . it is interesting to compute the differential conduc- 
tance of our DBMTJs, for the P majority channel. These 
are shown in figure [5l The differential conductance is de- 
fined as G = dl/dV, and in this work it is simply cal- 
culated numerically from the I-V curves. From Figure [5] 
it is seen that both DBMTJs show a conductance peak 
around V — 0.1 Volt. There is also a second conductance 
peak occurring at around V — 0.35 Volt for m=2 MLs 
and at around V = 0.5 Volt for m=4 MLs. Importantly 
the bias voltages at which these peaks occur are rather 
close to the experimental ones, even if our DBMTJs have 
significantly thinner in-between Fe slabs. 

A second intriguing feature of the DBMTJs is that the 
P majority current does not increase linearly with bias, 
in contrast to what happens for both the P minority and 
AP currents. As a consequence the TMR decay is not 
as severe as in the case of single-barrier junctions. In 
addition to that, it is important to report that the AP 
current for DBMTJs is significantly reduced with respect 
to that of the SBMTJs. For the simple SBMTJ inves- 
tigated here, this is always larger than the P minority 
current, except at very low bias where the resonance of 
Fe minority surface state plays the dominant role Q . In 
contrast, for double-barrier junctions the AP and the P 
minority currents are almost equal to each other, even at 
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FIG. 5: Differential conductance associated to the majority 
current in the P configuration for DBMTJs with m—2 MLs, 
m=4 MLs, as a function of bias. 



relatively high voltages. This feature increases the bias 
voltage at which the TMR starts to decay pushing Vi/ 2 
to higher biases. That is to say, our DBMTJs can sus- 
tain large TMR values even under the application of high 
bias voltages. As it was mentioned in the introduction, 
the reduction of the AP current is an indication of the 
SFE originated by the insertion of a magnetic slab. Wo 
then conclude that the spin-filter effect plays an impor- 
tant role in the TMR decay with bias in double-barrier 
junctions. 

By comparing the DBMTJs with the SBMTJ, it is 
clear that the TMR of DBMTJs does not only decay 
slower but it is also considerably higher than that of the 
SBMTJ with the same barrier thickness. This is due 
to the reduction of AP current in DBMTJs, originated 
from the spin- filter effect jl4| . Therefore, the SFE in 
DBMTJs increases both V 1/2 and the TMR. These two 
results, together with the TMR oscillations at low bias 
mentioned above, are the most important features of our 
DBMTJs. It would be interesting to study the I-V curves 



of DBMTJs with a very thick in-between Fe slab, since, 
according to our picture based on the SFE, the TMR 
decay with bias would be significantly reduced and the 
TMR values would increase. 



V. CONCLUSIONS 

In summary, we have compared first principles cal- 
culations of the I-V characteristics of single and dou- 
ble Fe/MgO (001) magnetic tunnel junctions in the 
phase-coherent spin-conserving transport regime. Our 
results are in semi-quantitative agreement with recent 
experiments and demonstrate that: (i) the TMR for 
double-barrier junctions is significantly higher than that 
of single-barrier, not only at low but also for finite 
bias; (ii) the TMR decay with bias voltage is slower 
in double-barrier than in single-barrier junctions; (hi) 
double-barrier junctions exhibit low bias TMR oscilla- 
tions. We have shown that the spin-filter effect is impor- 
tant in understanding the origin of the features (i) and 
(ii) , while the third aspect is related to resonant electron 
tunneling through a majority spin QWS formed in the 
in-between Fe slab. In our opinion the first two aspects 
are the most relevant for real DBMTJs, since they are 
rather robust with respect to changes in the junction ge- 
ometry. Thus, a possible route to obtain a very large 
TMR, that does not depend strongly on the bias volt- 
age, is that of fabricating double-barrier junctions with 
very thick in-between magnetic slabs. This is a poten- 
tially useful strategy for constructing spin-transfer ele- 
ment, where one needs to increase the current density 
and preserve a high degree of spin-polarization. 
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